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Abstract 



In this paper we generalize the Jaynes-Cummings Hamiltonian by making use of 
some operators based on Lie algebras su(l,l) and su(2), and study a mathematical 
structure of Rabi floppings of these models in the strong coupling regime. We show 
that Rabi frequencies are given by matrix elements of generalized coherent operators 



( quant-ph / 0202081 ) under the rotating-wave approximation. 

In the first half we make a general review of coherent operators and generalized 
coherent ones based on Lie algebras su(l,l) and su(2). In the latter half we carry out 



a detailed examination of Frasca ( |quant-ph /0 1 1 1 134 ) and generalize his method, 



and moreover present some related problems. 

We also apply our results to the construction of controlled unitary gates in Quan- 
tum Computation. Lastly we make a brief comment on application to Holonomic 
Quantum Computation. 



*E-mail address : fujii@yokohama-cu.ac.jp 



1 Introduction 



Coherent states or generalized coherent states play an important role in quantum physics, 
in particular, quantum optics, see Q and @. They also play an important one in math- 
ematical physics. See the textbook ||. For example, they are very useful in performing 
stationary phase approximations to path integral, [[|, 0, |J. 

Coherent operators which produce coherent states are very useful because they are 
unitary and easy to handle. The basic reason is probably that they are subject to the 
elementary Baker- Campbell-Hausdorff (BCH) formula. Many basic properties of them 
are well-known, see [|J or ||. 

Generalized coherent operators which produce generalized coherent states are also 
useful. But they are not so easy to handle in spite of having the disentangling one 
corresponding to the elementary BCH formula. In |7| and |14[ the author determined all 



matrix elements of generalized coherent operators based on Lie algebras su(l,l) and su(2). 
They are interesting by themselves, but moreover have a very interesting application. 
In [[[^] Frasca dealt with the Jaynes-Cummings model which describes a two-level 



atom interacting with a single radiation mode (see for a general review) in the strong 
coupling regime (not weak coupling one !) and showed that Rabi frequencies are obtained 
by matrix elements of coherent operator under the rotating-wave approximation. His aim 



was to explain the recent experimental finding on Josephson junctions JIT . 

This is an interesting result and moreover his method can be widely generalized. See 
also []n| for an another example dealt with in the strong coupling regime. 



In this paper we generalize the Jaynes-Cummings Hamiltonian by making use of some 
operators based on Lie algebras su(l,l) and su(2), and study a mathematical structure of 
Rabi floppings of these extended models in the strong coupling regime. 

We show that (generalized) Rabi frequencies are also given by matrix elements of 
generalized coherent operators under the rotating-wave approximation. We believe that 
the results will give a new aspect to Quantum Optics or Mathematical Physics. 

We also apply our results to the construction of controlled unitary gates in Quantum 



Computation in the last section. 

Lastly we discuss an application to Holonomic Quantum Computation, but our dis- 
cussion is not complete. 

2 Coherent and Generalized Coherent Operators 

2.1 Coherent Operator 

Let a(a') be the annihilation (creation) operator of the harmonic oscillator. If we set 
N = a^a (: number operator), then 

[N,at] = a) , [N,a\ = -a , [a\a\ = -1 . (1) 

Let H be a Fock space generated by a and a* , and {\n)\ n e N U {0}} be its basis. The 
actions of a and a) on TL are given by 

a\n) = y/n\n — 1) , a'\n) = \Jn + \\n + 1) , N\n) = n\n) (2) 

where |0) is a normalized vacuum (a|0) = and (0 1 0) = 1). From (0) state \n) for n > 1 
are given by 

|n) = ^tIO) . (3) 
These states satisfy the orthogonality and completeness conditions 

oo 

(m\n) = S mn , $3|n)(n| = l. ( 4 ) 

ra=0 

Definition We call a state 

|z) =e** f -**|0) =U(z)\Q) for 2GC (5) 

the coherent state. 

2.2 Generalized Coherent Operator Based on su(l, 1) 
Let us state generalized coherent operators and states based on su(l, 1). 
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We consider a spin K (> 0) representation of su(l, 1) C si (2, C) and set its generators 
{K + ,K.,K 3 } ((K + y = K_), 

[K 3 , K + ] = K + , [K 3 , K.\ = -K-, [K+, K_] = -2K 3 . (6) 

We note that this (unitary) representation is necessarily infinite dimensional. The Fock 
space on which {K + , K_, K 3 } act is TCk = {\K, n) \n e N U {0}} and whose actions are 



K+\K,n) = y/(n + l)(2K + n)\K,n + l), K_\K,n) = yfn(2K + n - 1)\K, n - 1), 
K 3 \K,n) = (K + n)\K,n), (7) 

where 1-^,0) is a normalized vacuum (K-\K,0) = and (K, 0\K, 0) = 1). We have 
written \K,Q) instead of |0) to emphasize the spin K representation, see ||. From (0), 
states \K,n) are given by 

\K,n)= \K,0), (8) 

/ n!(2iT) n 



where (a) n is the Pochammer's notation (a) n = a(a + 1) • • - (a + n — 1). These states 
satisfy the orthogonality and completeness conditions 

oo 

(K,m\K,n) = 5 mn , J2\ K > n )( K > n \ =1 «- ( 9 ) 

n=0 

Now let us consider a generalized version of coherent states : 
Definition We call a state 

\z) = V{z)\K,0) = e zK+ -~ zK ~\K,0) for z G C. (10) 

the generalized coherent state (or the coherent state of Perelomov's type based on su(l, 1) 
in our terminology). 

Here let us construct an example of this representation. First we set 

K + = \ (a^ , = \c? , K 3 = \ {a} a + i) , (11) 

then it is easy to check that these satisfy the commutation relations @. That is, the set 
{K + , K_, K 3 } gives a unitary representation of su(l, l)with spin K = 1/4 and 3/4, ||. 
Now we also call an operator 

S( z ) = e \{^) 2 -^} f or z e C (12) 



the squeezed operator, see the book ||. 



2.3 Generalized Coherent Operator Based on su(2) 

Let us state generalized coherent operators and states based on su(2). 

We consider a spin J (> 0) representation of su{2) C s/(2,C) and set its generators 
{J + ,J-,J 3 } ((J + )t = J_), 

[J 3 ,J+] = J + , [J 3 ,J-] = -J_, [J+,J_] = 2J 3 . (13) 

We note that this (unitary) representation is necessarily finite dimensional. The Fock 
space on which {J + , J_, J 3 } act is = {|J, n)\0 <n< 2J} and whose actions are 



J + \J,n) = y/(n + i)(2J-n)\J,n + l), J-\J,n) = s]n{2J - n + l)\J,n - 1), 
J 3 \J,n) = (-J + n)\J,n), (14) 

where \J, 0) is a normalized vacuum (J_|J, 0) = and (J, 0|J, 0) = 1). We have written 
| J, 0) instead of |0) to emphasize the spin J representation, see M. From flT4]), states 
| J, n) are given by 

\J,n) = ^=\J 1 0). (15) 
These states satisfy the orthogonality and completeness conditions 

2J 

(J,m\J,n) = 5 mn , ^2\J,n)(J,n\ = lj. (16) 

n=0 

Now let us consider a generalized version of coherent states : 
Definition We call a state 

\z) = W{z)\J,0) =e zJ+ -~ zJ -\J,0) for z G C. (17) 

the generalized coherent state (or the coherent state of Perelomov's type based on su(2) 
in our terminology). 

A comment is in order. We can construct the spin K and J representations by making 
use of Schwinger's boson method. But we don't repeat here, see for example @. 



3 Matrix Elements of Coherent and Generalized Co- 
herent Operators • • • [jlHj| 



3.1 Matrix Elements of Coherent Operator 

We list matrix elements of coherent operators U(z). 

The Matrix Elements The matrix elements of U(z) are : 



(i) n<m (n\U(z)\m) =e~-^ \ —(-z) m - n Lj m - n \\z\ 2 ), (18) 

V Tn\ 

(ii) n>m (n\U(z)\m) =e-^ zl \^z n - m Lj n - m \\z\ 2 ), (19) 



where L„ ^ is the associated Laguerre's polynomial defined by 



k / 
3=0 \ 



k + a\ xi 
k-j) j\ 



(20) 



In particular = is the usual Laguerre's polynomial and these are related to 

diagonal elements of U(z). 

3.2 Matrix Elements of Coherent Operator Based on su(l, 1) 

We list matrix elements of V(z) coherent operators based on su(l, 1). In this case it is 

always 2K > 1 (2K = 1 under some regularization). 

The Matrix Elements The matrix elements of V(z) are : 



n\m\ , 2-, k n+m 



(ii) n>m (K,n\V(z)\K,m) = J — - -—- K ^(l + NT K ~~ x 



where 



ks smh^ = _ (23) 
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The author doesn't know whether or not the right hand sides of (^T|) and (^) could be 
written by making use of some special functions such as generalized Laguerre's functions 
in (|20|). Therefore we set temporarily 

F {n ~ m) (x ■ 2K) = Vf-ir' j T(2K + m + n-j) x \i x ™-i (24) 

m { ] H r(2K)(m-j)\(n-j)\j\ [L + X)X [M) 

and Fj°\x;2K) = F m (x;2K). 

3.3 Matrix Elements of Coherent Operator Based on su(2) 

We list matrix elements of W(z) coherent operators based on su(2). In this case it is 
always 2 J e N. 

Matrix Elements The matrix elements of W(z) are : 



(i) n<m (J,n\W(z)\J,m) = J ^ (-R) m - n (l - \k\ 2 ) j ~^ x 

V 2jPn2jPm 

E*(-ir J (2J _ - _ - - Mm _ Mn _ m (i - ^\ 2 y^\ 2 r\ (25) 



(ii) n>m (J,n\W(z)\J,m) = \\ l'™' K n ' m {\ - \k\ 2 ) j -^ x 



2J r n2,J r rn 



where 



E-C-^Tor — ^ — Tj _ y( 1 - N 2 HN 2 r-M26) 

^0 (2J - m - n + j)!(m - - 



sin(lzl) ,. , n „ s 

k = , u z = cosfhz W (27) 



Here J2 * means a summation over j satisfying 2J — m — n+j > 0. 

The author doesn't know whether or not the right hand sides of (^) and (^Bj) could 
be written by making use of some special functions. We set temporarily 

m (2T\\ 

Fj n ~ m \x :2J) = Y: *(-l) m ^'7^ rV^ W( " (28) 

j=0 (2J - m - n + j)\(m - j)\(n - 

and Fj°\x; 2 J) = F m (x;2J). 
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4 Jaynes— Cummings Models in the Strong Coupling 
Regime 

In [|l^] Frasca treated the Jaynes-Cummings model and developped some method to 
calculate Rabi frequencies in the strong coupling regime. We in this section generalize the 
model and method, and show that Rabi frequencies in our extended model are given by 
matrix elements of generalized coherent operators under the rotating-wave approximation. 
This gives a unified approach to them. 

Let {o"i, a 2 , 03} be Pauli matrices and 1 2 a unit matrix : 

/ 



1 

1 



C"2 




03 




1 
1 



(29) 



The Hamiltonian adopted in [JT2 
(N) 



is 



A 



H N = ul 2 ® a) a + —a 3 <g> 1 + ga\ <g) (a 1 " + a) 



(30) 



where u is the frequency of the radiation mode, A the separation between the two levels 
of the atom, g the coupling between the radiation field and the atom. 
Moreover we want to treat the following Hamiltonians (our extension) 



(K) H K = ul 2 ®K 3 + —a 3 ®l K + ga 1 ®{K + + K_] 
(J) Hj = ul 2 ® J 3 + y03 ® 1j + g<ri ® (J+ + J-). 



(31) 
(32) 



To treat these three cases at the same time we set 

(N) {a\a,N} 
{L+,L_,L 3 } = | (K) {K + ,K_,K 3 } 
(J) {J + ,J-,J 3 } 

and 

H = H + V = lu1 2 ® L 3 + ^0-3 ® 1 L + g<J\ ® (L + + L. 



(33) 



(34) 



where we have written H instead of Hl for simplicity. 



Mysteriously enough we cannot solve these simple models completely (maybe non-integrable), 
nevertheless we have found these models have a very rich structure. 
For these (non-integrable) models we usually have two perturbation approaches : 
Weak Coupling Regime (0 < g < A) 



H = tol 2 <8> L 3 + — cr 3 <g> 1 L , V = go x <g> (L + + L_). 



Strong Coupling Regime (0 < A <C g) 

H = ul 2 <g> L 3 + go\ <g> (L + + L_), 



(35) 



(36) 



In the following we consider only the strong coupling regime (see for the weak one) 
First let us solve H which is a relatively easy task. 
Let W be a Walsh-Hadamard matrix 



111 
w = -7= 



then we can diagonalize cri by using this if as uj = \Va3W 1 . The eigenvalues of cti is 
{1,-1} with eigenvectors 



We note that 



|l)(l|=o 



-l)(-ll 



1 1 
1 1 
/ 

' v - 



w 



V 



1 1 



|A)(A| 



1 A 
\X 1 



W 



w 



I 1+A 
2 



w~ x . 



1-A 
2 



Then we have 



#0 = (W 7 ® 1l) (wl 2 ® L 3 + gcrs ® (£+ + £-)) 1 ® 1l 



, ujL 3 + q(L + + L_) . 

(W®1 L ) | | (w- 1 ®^). 

uL 3 -g(L + + L_ 
= |1)(1| <g> {cuL 3 + #(L+ + + | — 1 > < — 1 1 <g> {wL 3 - g{L + + 

= ^|A}(A|®{a;L 3 + A^(L + + L_)} 

A 

= E l A X A l ® {e-^( L +- L -) (QL 3 )e^- L -^} 
x 

= E (l A ) ® e-^^-)) (OLb) ((A| ® e*^- £ -)) (37) 

A 

where we have used the following 
Key Formulas For A = ±1 we have 

(iV) uoa)a + A(?(a f + a) = fie-^ (at " a) ^JV - ^ e ^ (at " a) 

where = u>, x = 2g/u>, (38) 
(K) uK 3 + Xg(K + + K_)=tt e-^ K +~ K -^K 3 e^ K +~ K ^ 



where Q = u^Jl - (2g/u) 2 , x = tanh -1 (2#/u;), (39) 
(J) u.J 3 + \g(J + + J_) = e - 2 ^ (J +- J - ) J 3 eT( J +~ J -) 



(Eigenvalues, Eigenvectors) 



where = w^/l + (2#/cj) 2 , ac = tan _1 (2#/u;). (40) 

The proof is not difficult, so we leave it to the readers. That is, we could diagonalize 
the Hamiltonian H Q . This is two-fold degenerate and its eigenvalues and eigenvectors are 
given respectively 

' (AO um-£, | A) ®e-iH at - a V) 
(K) n(K + n) 7 \\}®e-^ K +- K -^\K,n) (41) 
(J) n(-J + n), |A) ®e-^( J +- J -)|J,n) 

for A = ±1 and n£NU {0}. For the latter convenience we set 

Eigenvalues = {E n }, Eigenvectors = {|{A, n})}. (42) 
Then (|37D can be written as 

^o = EE^n|{A,n})({A,n}|. (43) 

A n 
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Next we would like to solve the following Schrodinger equation : 

z-^tf = H^= (h + -cr 3 1 L ) 



(44) 



where we have set ft = 1 for simplicity. To solve this equation we appeal to the method 
of constant variation. First let us solve 



d 
at 



which general solution is given by 



(45) 



*(t) = E7 (t)*o = e- itH °y 



(46) 



where \l/o is a constant state. It is easy to see from ([43] 



u (t) = e- w ° =5;E e ^l{^}>({^}l- 



A n 



We have only to solve this equation. If we set 



H F =^U \<r 3 ®l L )U , 



then we have easily from 



A 



(47) 



The method of constant variation goes as follows. Changing like — ► $?o(t), we 
insert fl46f) into (^). After some algebra we obtain 

.d . A 



Hf = T E E e^-^({A, m}|(a 3 ® U)|{/i, n}> |{A, m}}({//, n}| 

A,/z m.,n 



itn(m-n) / / | e Aa;(i + -I/_) 



n» |{A,m})({-A,n}| 



where we have used the relation (A|o"3 = (— A|. Remind that |n)) is respectively 



(N) |n) 
(K) \K,n) 
(J) |J,n) 



(48) 



(49) 



(50) 
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In this stage we meet matrix elements of the coherent and generalized coherent 
operators e Xx ( L +- L -) j n section 3 (z = z = Ax). 
Here we divide Hp into two parts 



where 



Hp = Hp + Ht, 



Hf =^EE<He Ax,(L+ - L -V>> \{X,n})({-X,n}\, (51) 

z A n 

Hp" = -EE e im( - m ' n \(m\e Xx{L+ - L - ) \n}} |{A, m})({-\, n}\. (52) 



Noting 



((nle^+^-V)) = ((nle^^+^-V)) 



by the results in section 3, ifp can be written as 

Hf = y E ((n|e^- L ->|n» {|{l,n}>({-l,n}| + |{-l,n»({l,n}|} , 
from which we can diagonalize if^ as 

^' = vEE((^ (i+ " L - ) N)^ IfolMXfolMI (53) 

Z n (j 

if we define a new basis 

\{a^ n }) = -L( ff |{l,n}> + |{-l,n}», a = ±1. 



These states can be seen as so-called Schrodinger cat states, |T3[ . From these we have 



I{l,n}> = ^{|{l^»}>-|{-1^»}>}, 
|{-I,n}) = -^{|{1,^}) + |{-1,^})}. 



Inserting these equations into ( p2j ) and taking some algebras we obtain 

V = 7 E E^ (m - n) {((^|e^ + - L -V))| lfoiU><{</^»}| + 



(( m \e-^- L -^\n)}j IfaMHVMll. (54) 
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For simplicity in we set in the following 



E n>a = ^a((n\e x ^- L -^\n)), 



(55) 



then 



En 



-ere 



(N) 

(K) ^(Sr( 1 + \K\ 2 )- K ~ n F n {\K\ 2 : 2K) where re = sinh(x) (56) 



(J) f - \K\ 2 ) J ~ n F n {\K\ 2 : 2J) where re = sin(x) 



from fl41|) and the results in sectin 3.1. Now let us solve 



For that if we set \P (^) as 



*o(t) = EE e "^"' CT ^ CT (t)|{o-,^}), 



(57) 



then we have a set of complicated equations with respect to {a nj(J }, see |12j. But it is 
almost impossible to solve them. Therefore we make a daring assumption : for m < n 



(5f 



That is, we consider only two terms with respect to {n\n > 0}. After some algebras we 
obtain 

ij t a m „ = ^e-^'-^^e^— ) {((m\e^-^\n))^ + ({m\e~^-^n)) a -\ a n ^ 

(59) 

But we cannot still solve the above equations exactly (see Appendix), so let us make 
so-called rotating-wave approximation. The resonance condition is 



{E na > - E m>a ) + {m - n)tt = 



E na > -E m>ff = (m-n)fi (60) 



for some a and a , and we reject the remaining term in (|59|) . Then we obtain simple 
equations : 
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Interband Transition Case (er 7^ a) E n> -, a — E my(7 = (m — n)Q 



.d 
at 



|{(He^- i -)|n))|-((m|e-^- L -)|n))|}a n , 
^er((m|sinh (x(L + - LJ)) |n))a„_ (T , 



* 7 1 ^n,—a 

at 



A 
2" 



£ {-((n|e^+~ L -V))| + ((^|e- x(at - a V))|} <W 



A 



cr((n|sinh — LJ)) |m))a OT)(J . 



Intraband Transition Case (er = er') i£ n)<7 — i? mi<7 = {m — n)Vl 



d 



i—a r 



dt 



— |((m|e 

^er((m|cosh(x(L + - LJ)) |n))a nj<7 



^ 1 -i-)| n ))| + (( m |e-^ + -^)| n ))|^ /((f7 



.d 

z— a. 
dt 



n,cr 



A 
2" 
A 
2" 



E|((n|e^- L -)|m))| + ((n|e-^- L -)|m))|}a m , CT 

0"' 

cr((n|cosh (x(L + - LJ)) |m))a TOi<J . 



then 



(61) 



(62) 



For simplicity we set 



7e = A((n|sinh(x(L+ -L_)) |m)), ft' = A((n|cosh (x(L+ - LJ)) \m)), (63) 



A((m|sinh (x{L + - LJ)) \n)) = -11, A((m|cosh (x(L + - LJ)) \n)) = Tl' . 

These are two Rabi frequencies as shown in the following. It is important that Rabi 
frequencies in our models are given by matrix elements of coherent and generalized 
coherent operators ! 

By making use of the results in section 3 and (0H), (|3"5|), (HO) we have 



(AT) 



* = f \/S (^)"" m e-^L m (— ) (£) {1 - (-1)—} 



ft' 



A / to! 
2 V n\ 



(64) 
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K 



n\m\ 



2 V (2K) n (2K), 
A 



-K 



1 + AC ) 2 



Fj n - m \\K\ 2 : 2K){1 - {-l) n - m } 



?. \ m. \ 



(J) 



2 V 2 J Pn 2 J Pm 



2jPn2jPm 



1 + \K\ 2 y K ~^Fj n - m \\K\ 2 : 2K){\ + (-I)™"" 1 } 

V W 

1 - \K\ 2 ) J -^Fj n - m) {\K\ 2 : 2J){1 - (-1)™—} 



(65) 



(2K) n (2K) m " 

where k = sinh fx) with x = tanh 



- H 2 ) J -^F m (n - m) (M 2 : 2J){1 + (-1)™-™} 



where k = sin (x) with x = tan 1 (^~^J 
From these we find a constraint between m and n : 

Interband Case n — m = 2N — 1 =^ n = m + 2N — 1 for JVgN, 



Intraband Case n — m = 2N 
Now let us solve (|6l| ) and (|62|). 



n = m + 2iV for iVGN. 



/ 



a„ 



( 



( 
\ 
( 



-o 



7v 



-af 



a 



'A 
J 



so their solutions are given by 

a m,a(f) \ [ COs(yt) iaSm(J^t) 



\ a n - a {t) 



\ an,a(t) 



icxsin(yt) cos(-jt) 



<w(0) 

\ On ,-ff(0) 



cos(^t) — iasm(^-t) 
-zcrsin(^t) cos(^t) 



V a„, CT (0) 



(66) 



(67) 



(6* 



We have obtained some solutions under the rotating-wave approximation. Now it may 



be suited to compare our results with a recent experimental finding in [11], but this is 



beyond our scope. See |12 



Let us conclude this section by a comment. Our ansatz (|58|) to solve the equation 
is too restrictive. We want to use (|57D to solve the equation, but it is very hard at this 
stage. 
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Problem Find more dynamic methods ! 



5 Quantum Computation 

Let us reconsider the results in the preceding section in the light of Quantum Computa- 
tion. Remind once more that the following arguments are based on the rotating-wave 
approximation. 

Interband Case {a = 1) 





( 




\ 


/ 











1Z 
2 


\ 


/ 




\ 


d 




-1 


























i— 

dt 




0"n,l 






















a n ,i 






\ 


tin -1 


/ 


V 


n 

2 











/ 


V 




J 



The solution is 



/ 


a m ,i(t) 


\ 


( 


cos(ft) 








isin(yt) 


\ 


/ 


a m ,i(0) 


\ 




a m -i(t) 









1 












a m ,-i(0) 






a n ,i(t) 












1 









a„,i(0) 




V 


a n ,-i(t) 


J 


V 


isin(-jt) 








cos(ft) 


/ 


V 


On-l(O) 


/ 



(70) 



Interband Case (a = — 1) 





( 


fl m,l 






( 














\ 




/ 






d 




"m,-l 














n 

2 











a rn-l 




I— 

dt 




On,l 











K 
2 














0>n,l 






{ 


0>n,-l 


J 




{ 














/ 






Q>n-1 


) 



(71) 
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The solution is 
/ 



a m ,i(t) 
S-i(') 

\ a n ,-i(t) j 



10 

cos(ft) -isin(fi) 
-isin(ft) cos(ff) 



\ 



' o m>1 (0) ' 

Om-l(O) 
On,l(0) 

V Q n-i(°) y 



Intraband Case (a = 1) 



d 
i — 
dt 



O-m-1 



\ 



1 0^0^ 




The solution is 
/ 



v a n _i y yooooyy a„_i y 



^ o 

2 w 







a m ,i(0 
a m -i(t) 
a„,i(t) 

i(0 







-isin(^-i) 











— isin(^-t) 




Intraband Case (er = — 1) 



cos(^- 




^t) 



1 



7 V 



Om,l(0) 

o m -i(0) 

On,l(°) 
On-l(O) 



\ 



/ 



; d 
'dt 



V a n,-l J 



The solution is 
/ 



a m ,i(t) 
a m -i(t) 
a n ,i(t) 





-^ 



0-^0 



0"m,-l 

\ a n -i ) 



10 

cos(^) isin(^-f) 

10 

^ isin(^-t) cos(^t) j 
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Om,l(0) 

a m _i(0) 
a„,i(0) 

V an-i(o) y 



If we can identify Q58D with an element in two-qubit space 



a m , 1 (*)|00) + a m ,_i(t)|01) + o n ,i(t)|10) + a n ,_i(0|ll) e C 2 ® C 2 



(77) 



where C 2 = Vect c {|0), |1)}, then the solutions ([7DJ) , (|7J), (0), ([7BJ) are kinds of con- 
trolled unitary operations (gates) which play a crucial role in Quantum Computation, 
see for example |fTB| . For example, ([7E|) is just one of controlled unitary gates expressed 
graphically as 



We note here that controlled unitary gates above are written down as 



C-Unitary 



for 



U 



^ 1 ^ 



Mn W12 

10 

y0 U 2 1 o U 22 J 



(7E 



till Uu 



U2\ U 2 2 



e U(2). 



A comment is in order. [0 (and (T^|) is considering the same subject. However the 
authors in them treated it in the weak coupling regime, while we treated it in the strong 
coupling regime. It is very interesting to investigate a deep relation (connection) between 
them. 
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6 Discussion 



One of motivations of this study is to apply our results to Holonomic Quantum Computa- 
tion developped by Italian group (Pachos, Rasetti and Zanardi) and the author, see . 
23, [0], @ and J2§— JZ7| and recent 



In this theory we usually use the effective Hamiltonian of a single-mode field of Kerr 
medium 

H Q = XN(N-1), N = a i a where X is a constant (79) 
as a background and the real Hamiltonian is in one-qubit case given by 

H(z, w) = W(z } w)H W- 1 (z, w) (80) 

where W is a product of coherent operator U(z) and squeezed one S(w) in section 2. In 
the above Hamiltonian Hq the zero-eigenvalue is two-fold degenerate whose eigenvectors 
are |0) and |1). We set \vac) = (|0), |1)). Then we can construct a connection form A on 
the parameter space {(z,w) G C 2 } as 

A=(vac\W~ 1 dW\vac) (81) 

from flSOl) where d = dz-^ + dw-J^. By making use of this connection we can construct 
a holonomy group Hol(A) (C U(2)) which is in this case equal to U(2). In Holonomic 
Quantum Computation we use this holonomy group as unitary operations in Quantum 
Computation. The point at issue is that we use not full property of the Hamiltonian but 
only property of the zero-eigenvalue. 



By the way, the Hamiltonian Hp in (19) 



H F = ^-Uo^ias® 1 L )U where L — (N) or (K) or (J) 

is very similar to fl8"U|). This system is always two-fold degenerate. Then a natural question 
arises : 

Problem Is it possible to perform a holonomic quantum computation by combining 
the systems {(N), (K), (J)} ? 

18 



This is a very interesting and challenging problem. 
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Appendix 



On Equations (|59|) 



Here let us write down full equations of fl5*9"|) with matrix equation form : 



'dt 



where 




for k = m, n 



and 



A = A(t) 



( 'BLjt(-E n , 1 +E mA +n( : m-n)) _K it{-E n - 1 +E m , 1 +Q.{m-n)) 

2 2 

TZ it(-E rl ^+E m - 1 +Q(m-n)) Tl' it(-E n .-i+E m .--,+n(m-n)) 

\ 2 C 2 C 

'ELpit(-E nA +E m ,i) _K p it(E nA +E m> i) 

jtn(m-n) I 2 c 2 c 

TZ it(-E„,i-B m>1 ) _BL p it(E nA -E m .i) 

2 2 



(82) 



(83) 



because _<t = —E^. 

We can give flg2p a formal solution by infinite series (called Dyson series in Theoretical 
Physics). Then we meet secular terms. 
For example let us consider the following simple equation : 



A a = e ^ a w ith a(0) = a 
dt 



The solution is given by 



a(t) 



c e 



exp (^) - ^ 
uo = 
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That is, we meet the secular term. 

By the way, we have known how to handle (simple) secular terms called Renormaliza- 



tion Group Method (Approach), see []30| for a general introduction. 

Frasca in |H] has applied this method to the above equation. The conclusion is 
interesting, but seems to be rather involved. We are now reconsidering his approach. 
Therefore let us present 

Problem Solve this matrix equation completely ! 
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